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Abstract 

We consider interactions of fermions with the domain wall bubbles produced 
' during the first order phase transitions. New exact solution of Dirac equations and 
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^1 During these last few years, a considerable amount of works has been 

devoted to the analysis of possibility of generation of baryon asymmetry 



of the universe on the electro- weak scale (see for example |], @, |3j]). In 
these models usually considered first order phase transitions, which can 
be described in terms of nucleation of bubbles of "true" vacuum (with an 
inner vacuum expectation value of Higgs field d / 0) appearing in the 
preexisting " false" vacuum (with v = throughout) and growing until the 
Universe is completely converted to the broken phase. So it is important 
to investigate interaction of particles from cosmological plasma with new 
phase bubbles. In this paper the new exact solution of Dirac equations is 
obtained for realistic bubble wall potential. 

The tree potential of the Higgs field has the form: 

v = ^V-A (i) 

where A is the self-coupling constant of the scalar field. Once the new 
phase bubble has grown to macroscopic size the curvature of its wall can 
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be neglected and can be considered as a flat. The world is divided in two 
zones, on the left-hand side of the wall, say, v = 0, while on the right-hand 
side — v 7^ 0, and, thus due to Yukawa coupling quarks and leptons acquire 
their masses. In this frame it is assumed that the bubble expands from 
right to left. The analytic solution of Higgs field equation for a bubble wall 
with its normal along the z-axis and position at z = o is 

' z" 
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(2) 



(4) 



where _ 

is the width of the wall and M# = y/Xv/2 is the Higgs mass. We want to 
emphasize that the domain wall profile (|2j) follows directly from the field 
equations for the Higgs field in difference to various "written by hand" non 
realistic wall-profiles which have been considered |2], [3f . 

For the simplicity let us consider only one fermion flavor. The fermion 
mass 

Tfl ( Z 

M(z)=g<p(z) = — [l + th\^ 

(where g is Yukawa coupling constant) is a function of position with the 
boundary conditions 

M(-oo) = 0; M(oo) = gv = m, (5) 

i.e. it is zero outside and nonzero inside the bubble. 

Let us now consider scattering problem. Assume that incident particle, 
coming from the left, in the domain wall frame interacts with the wall. 
This puzzle is very much like the quantum mechanical scattering of Dirac 
particles from a potential barrier The difference from the quantum 
mechanical case is that an exterior field appears in the mass of scattered 
particle. Thus we are solving Dirac equation with varying mass. In order 
to extract the reflection coefficients we need to find plane wave solutions 
corresponding to incoming and outgoing particles moving along z-axis. 
The components of the momentum parallel to the wall are conserved in 
the scattering process and are unimportant in a qualitative discussions. In 
order to solve the problem of reflection from the wall one has to specify 
the boundary conditions. Particles moving in the unbroken phase have a 
wave function proportional to 
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where the upper sign corresponds to incident particle moving towards the 
wall, while the lower sign - - to the reflected particle. The wave function 
of particles transmitted through the wall is proportional to 

Here p = \/E 2 — m 2 , q = \E\. The boundary conditions for antiparticles 
are the similar as for particles. 

We shall work in chiral basis (the notation of 7-matrixes the same as 
in Ref. || ). For the particles with positive energy Dirac equation with a 
position-dependent mass (H) can be written in components in the form 

1— ^1 + Eik + M(z)^ 3 = 0, 

i _ E fo _ M{z)il)i = 0, (8) 
az 

where ipi and ips correspond to the left- and the right- handed particles. 
Introducing a new variable 

and searching the solution of (|8]) in the form 

^ m ^y- iSp/2 (i-y) ±iSE/2 v m (y), (10) 

(where the upper sign in the second exponent corresponds to ipi, while the 
lower one — to ^3) one yields the equation: 

{y(l - y)j^ + [(1 - i6p) - (1 - iSp ± iSE)y]j- - (11) 

-^5 2 E(E T p)}r]i(3) = 0. 

(We want to note that in the first exponent in (pD|) we could choose the 
opposite sign, however in that case we would yield the outcoming plane 
wave at +00 with incorrect sign of momentum p) . Equation ([T2]) is the one 
for Hypergeometric functions M and has a solutions 



Vi(3) = F 
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-(m — ip± iE), — (— m — ip ± ie), 1 — iSp, y 



(12) 



Now let us consider an asymptotical behavior of functions i/ji and ip%. In 
the area z —> —00, i.e. when y — ► 1 and (l — y)—> e 2z ^ s , using the features 
of hypergeometric functions we obtain: 

^1(3) - ^1(3)^- (13) 
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Thus ipi corresponds to the incident (coming from the left) left handed 
massless particle with momentum equal to \E\ while ^3 - - to the reflected 
particle of the right chirality. 

Constants ^.1(3) are expressed by Gamma functions in the following way: 



A(3) 



r(i - iS P ) ■ r(i ± iSE) 



r M-m - ip =f iE) + 1] -V 



ip =p iE) + 1 



(14) 



In the broken phase, when z —> 00, i.e. when y — > e 2z > 6 , (1 — y) —> 1, 
we have only the wave corresponding to transmitted massive particle: 



^1(3) 



(15) 



with momentum p = \J E 2 — m 2 . 

To find the reflection coefficient we must compute the ratio of reflected 
flux of particles to the incident one. Substituting into expression for 
flux we yield: 



J z = z^tV^ = VWi - = \A Y \ 2 - \A 3 \ 2 = J\ 
The reflection coemcient is: 
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Substituting expression ( |T^) in this formula and using the features of 
Gamma functions 



r(i + x + %Y) ■ r(i + x - iY) ■ r(i - x + zT) • r(i - x - %y) = (is) 

2tt 2 {X 2 + Y 2 " 



we have 



sh{2irY) - cos(2ttX) 
E + p sh[d7v(E — p)} — cos{5 , Km) 



^ _ J-J ^ p ^ onyunyj^ jj j j ^uoyvnni) 

E — p sh[5ir(E + p)] — cos{5nm) 

Thus in this letter we derived an analitical expression of the reflection 
coemcient. This expression can be used, for example, in calculations of 
baryon asymmetry of the universe. 
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